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Abstract 

This paper covers the propulsion system component modeling and controls development 
of an integrated nonlinear dynamic simulation for an inlet and engine that can be used 
for an overall vehicle Aero-Propulso-Servo-Elastic (APSE) model. The focus here is on 
developing a methodology for the propulsion model integration, which allows for controls 
design that prevents inlet instabilities and minimizes the thrust oscillation experienced by 
the vehicle. Limiting thrust oscillations will be critical to avoid exciting vehicle aeroelastic 
modes. Model development includes both inlet normal shock position control and engine 
rotor speed control for a potential supersonic commercial transport. A loop shaping control 
design process is used that has previously been developed for the engine and verified on 
linear models, while a simpler approach is used for the inlet control design. Verification of 
the modeling approach is conducted by simulating a two-dimensional bifurcated inlet and a 
representative J-85 jet engine previously used in a NASA supersonics project. Preliminary 
results are presented for the current supersonics project concept variable cycle turbofan 
engine design. 


N omenclat ur e 


2DB 

Two-Dimensional Bifurcated inlet 

HSR 

High Speed Research program 

A 

Cross-sectional area 

K 

Experimentally obtained bleed flow constant 

AIP 

Aerodynamic interface plane 

L 

Atmospheric disturbance wave length 

APSE 

Aero-propulso-servo-elasticity 

LAPIN 

Large Perturbation Inlet simulation 

ASE 

Aero-servo-elasticity 

M 

Mach number 

C inlet 

Inlet controller transfer function 

N+3 

Three generations forward from current state of the art 

CFD 

Computational fluid dynamics 

NPSS 

Numerical propulsion system simulation 

F 

Computational fluid dynamics flux vector 

P 

Pressure 



Q source 

Heat addition 
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Q 1 D Quasi-one-dimensional 

QFT Quantitative feedback theory 

R Gas constant 

S Inlet source vector 

S' Inlet bleed source vector 

T Temperature 

V Volume 

VCE Variable cycle engine 

W Inlet state vector 

e CFD internal energy term 

m Mass flow rate 

8 Laplace operator 

t Time 

u Velocity 

x Length 

A Change in variable 

a Artificial viscosity term 


7 Ratio of specific heat 

p Density 

Subscript 

1 Inlet plenum station 1 

2 Inlet plenum station 2 

c Corrected flow parameter 

n Station location 

s Static flow condition 

t Total flow condition 

v Engine volume gas parameter 

P Plenum 

Superscripts 

engine Location of the first engine volume 

inlet ex u Location of the aerodynamic interface plane 


I. Introduction 

U nder the NASA Fundamental Aeronautics Program, the Supersonics Project aims to overcome the 
obstacles associated with supersonic commercial flight. Current work on this project is to develop the 
technologies to allow for practical supersonic commercial flight at approximately Mach 2. The proposed 
vehicles are long, slim body aircraft with the potential for pronounced structural vibrations known as aero- 
servo-elastic modes (ASE ). 1,2 When coupled with propulsion system dynamics, the modes are known as 
aero-propulso-servo-elasticity (APSE). The APSE considerations can lead to design challenges pertaining to 
aircraft performance such as aircraft ride quality and stability. Furthermore, other disturbances upstream 
of the inlet generated by atmospheric wind gusts may also affect ride quality and stability. To study the 
propulsion component of APSE, an integrated nonlinear model was developed that includes both the inlet 
and engine propulsion system dynamics. 

The modeling approach used for the component dynamic models has previously been developed 3,4 and 
verification of inlet and engine integration modeling approach leverages work conducted under the NASA 
High Speed Research (HSR) project for an integrated mixed compression inlet and J-85 turbojet engine . 5 
The previous integrated inlet and engine was developed under HSR in the FORTRAN computer language. 
A legacy inlet code, LArge Perturbation INlet (LAPIN), and engine code the Aerodynamic Turbine Engine 
Code (ATEC ) 6 ’, 7 were used for the model integration. LAPIN is used as a key verification tool in the 
current development of the inlet model. Other integrated propulsion system simulations have been created 
by Gamble, Numbers, and Giannola 8,9,10 . The distinction in the work presented here is that the main goal 
is to provide a platform for control algorithm testing, in the context of changing design concepts for the 
supersonic propulsion system and integration into an ASE vehicle model. Under this goal the objective was 
to have a highly modular code built on a platform that is familiar to controls engineers. To this end the 
MATLAB® Simulink® software was chosen. 

The purpose of the supersonic inlet is to supply the engine with the required mass flow at a subsonic 
speed with high total pressure and minimum distortion. The inlet is modeled using a quasi-one dimensional 
(Q1D) implementation of the Euler equations. The engine comprises another major component of the overall 
propulsion system and is connected just downstream of the inlet flow path. Currently, the engine model uses 
a single lumped volume for each of the major components, such as the compressor, combustor, turbine, and 
afterburner-nozzle. The propulsion system controller is a third major component to the simulation. To study 
propulsion system controller effectiveness in maintaining the desired ride quality and stability, atmospheric 
turbulence models 11,12 are used to investigate the impact of free stream flow disturbances on the propulsion 
system. 

The engine control design approach is based on classical loop shaping ; 13 and quantitative feedback theory 
(QFT) methodologies laid out by Houpis . 14 The latter accounts for design specifications in terms of bounds 
in a loop shaping development approach . 15 The loop shaping approach used here is a linear controls design 
methodology that requires the nonlinear propulsion models to be linearized prior to the development of 
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control algorithms. The inlet control algorithm development used a Proportional Integral Derivative (PID) 
design that will be discussed in more detail subsequently. The developed algorithms are later tuned for 
implementation in the nonlinear model. 

The overall paper objectives are as follows: 

• Provide an overview of a propulsion system model for a supersonic vehicle, with the focus on integration 
of the inlet and engine. 

• Illustrate the key features of the propulsion system model to support developing and designing control 
algorithms. 

• Provide an example control case of using an inlet bypass door and an engine fuel injector to dampen 
out atmospheric disturbances, thus limiting thrust oscillations and allowing for stable operation. 


II. Model Overview 

The overall APSE simulation block diagram can be seen in Fig. 1. The propulsion system of the overall 
model, highlighted in Fig. 1, is comprised of the nonlinear inlet, engine, and propulsion controllers. ASE 
vehicle impacts are neglected in this study, but are illustrated by dashed lines in Fig. 1 to provide an overview 
of how the propulsion system fits into the larger APSE simulation architecture. All of the propulsion elements 
are modeled using unsteady conservation equations, described in the following sections. This allows for the 
investigation of inlet stability and thrust oscillations. To run the simulation, free stream static conditions 
of pressure, temperature, and Mach number are applied to the inlet model. The inlet model then calculates 
an exit total pressure and temperature, while controlling the terminal shock to maintain stability. Airflow 
properties at the inlet exit or aerodynamic interface plane (AIP) are used as the engine face input conditions, 
which establishes the key propulsion system’s integration point that is required for the two models. The 
engine then uses the input conditions, its performance maps, and a modified version of the Euler conservation 
equations to calculate the engine output thrust. The integrated model allows the engine simulation to initiate 
an AIP disturbance that impacts inlet shock stability. 

The overall goal of the ASE project 
is to integrate a high fidelity vehicle and 
propulsion system model to investigate 
aero-elastic effects. The work presented 
here is a step along that path, provid- 
ing a dynamic propulsion system simula- 
tion that is suitable for controls devel- 
opment and testing. Previous work in 
the area of APSE has focused on either 
the vehicle or the propulsion system, but 
rarely both, unless considering a hyper- 
sonic vehicle where the vehicle itself pro- 
vides compression of the flow feeding into 
the propulsion system. 16 18 



Figure 1. Block diagram of overall APSE system model with 
propulsion system highlighted. 


III. Analysis 

A. Inlet Model 

The inlet modeling methodology verification in this study was done using a simulation of the HSR 
project’s two-dimensional bifurcated (2DB) mixed compression inlet. 19 Work is currently ongoing to tran- 
sition the current design to an external compression inlet. However, the mixed compression design will be 
used for the remainder of this study to illustrate the modeling approach. 20 A mixed compression inlet is 
more complex due to the need to control the normal shock downstream of the inlet throat. If the normal 
shock moves upstream of the throat, it will be ejected from the inlet. This is known as an inlet unstart 
and results in a sudden loss of thrust. A schematic of the 2DB inlet can be seen in Fig. 2, which illustrates 
the desired region to control the normal shock and the bypass door actuator that controls the normal shock 
position by altering the mass flow. 

The inlet model developed using MATLAB Simulink software uses computational fluid dynamics (CFD) 
based numerical techniques to integrate the Q1D form of the Euler equations to model the dynamic flow 
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properties of the inlet. A schematic of the modeling approach can be seen in Fig. 3, where over a small 
distance, dx, the fluid properties change by a differential amount. The model assumes a calorically perfect 
gas. In addition, the code includes simple models of mass flow bleed and bypass. The bypass door is modeled 
as a choked bleed, since the main purpose of the bypass door is to release mass from the inlet. The bypass 
and bleed models will remove mass from the inlet core flow as illustrated in Fig. 3. 


Point of Low Risk, 
Low Performance 



Figure 2. 


2DB inlet schematic illustrating shock patterns and the location of the bypass door near the AIP. 


A brief synopsis of the modeling approach is provided here, 
without consideration for bleed or bypass as these will be dis- 
cussed in the following section. The approach uses the conser- 
vation equations of continuity, momentum, and energy, Eqs. 
(1) to (3), respectively: 21 


d(puA) 

dt 

dt 


9(pA) , djjmA) = 
dt dx 

d[A(pu 2 + P)} _ d A 
dx dx 


d[A(^ + ^-)] 

+ dx 0 


(1) 

(2) 

( 3 ) 


dx 


A 

u 

P 

P 



A+dA 

u+du 

P+dP 

p+dp 


Figure 3. Q1D flow representation of differ- 
ential volume. 


The above governing equations are implemented in the simulation assuming that ^ = 0. This requires 
the inlet geometry to be fixed during the simulation. This is acceptable because geometry conditions change 
very slowly when compared to supersonic flow. The resulting state, flux, and source vectors are provided 
below, Eqs. (4) to (6): 


P 


W = 


pu 

P I pu 2 
_7-l 2 _ 


( 4 ) 


F = 


pu 

P + pu 2 

Puj , pu 3 
_7-l 2 _ 


( 5 ) 


0 


s = 



The overall governing equation is expressed in vector form in Eq. (7): 


(6) 
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d W 

dt 


1 dAF S 
A dx A 


( 7 ) 


In the model, Eq. (7) is integrated numerically using a Runge-Kutta time marching scheme with central 
differences of the spatial terms. To increase stability in the vicinity of shocks, an artificial viscosity term is 
included to dampen out non-physical oscillations. This is used to smooth out numerical oscillations about 
the discontinuity of the normal shock. 21 The coefficient, a, is carefully chosen such that dampening occurs 
but the associated diffusion does not degrade the formal accuracy of the numerical scheme. The numerical 
viscosity term is essentially the equivalent of adding diffusion to the governing equations as shown in Eq. (8) 


d 2 AW 


a 


dx 2 


The governing equation is then updated to become Eq. (9): 


(8) 


c)W 1 dAF-a^jW s 
dt A dx + A 


(9) 


Bleed and Bypass Model 


Bleed zones and bypass gates are implemented throughout the inlet to control and maintain the position of 
the normal shock. The bleed sections are passive and add stability for the normal shock, while the bypass 
door is placed near the AIP for closed loop control of the normal shock position. Bleed zone gas states are 
calculated using the total pressure and total temperature of the grid point with an associated plenum that 
can fill to reduce the mass flowing in the main inlet duct section as illustrated in Fig. 4. The bleed zone 
plenums are required since the bleed flow from the main duct is not choked (Station 1). This allows flow to 
move in and out of the main duct from the plenum. The exit of the plenum into the free stream is assumed 
to be choked (Station 2). Each bleed section has its own plenum. The mass accumulates in the plenum, 
where the flow velocity is assumed to be zero. 22 

The governing equations used for the bleed plenum models 
are continuity Eq. (10) and energy Eq. (11), which are tied 
together by the state equation, Eq. (12). These equations are 
used to solve the gas properties at the stations referred to in 
Fig. 4. 22 




^ = TTTirMi - i - (7 - 1 )T P fn 2 ] 


Vppp 


Pp = ppRTp 


( 10 ) 

( 11 ) 

( 12 ) 



Figure 4. Bleed model schematic with as- 
sociated plenum volume, where the flow is 
unchoked at station one and choked at sta- 
tion two. 


Bleed flow rates are determined using an experimentally found variable discharge coefficient K. Equation 
(13) is used to determine the flow rates at station 1 and station 2, along with the gas parameters at the 
given station. 22 


fn = K 1 

ViRTt [i + 3±i ] 5(^i) 

The variable K is a function of bleed hole area, boundary layer, and Mach number. At station 1 the discharge 
coefficient is also a function of the pressure ratio across the bleed hole to account for the unchoked mass 
flow. Alternatively, an unchoked mass flow equation could be used at station 1, such as that found in Hill. 23 

The bypass door is modeled such that the flow out of the bypass door from station 1 is assumed to be 
choked. This eliminates the need for a bypass door plenum model. Additional source terms due to the bleed 
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and bypass are illustrated in Eq. (14). The new source is applied to the overall inlet governing equations by 
adding to the previous source vector Eq.(6). 22 


S' = - 


nii 

A n Ax n 


niiu n 
A n Ax n 


nii 


\ -7 Prt 
(7 -l)pn 


A n Ax n 


(14) 


Inlet Model Verification 


To verify the dynamic behavior associated with the present Simulink inlet simulation, the LAPIN flow code 
is used to represent the “true” solution. 22 LAPIN is a sophisticated high-resolution, Q1D code capable 
of modeling core flow instabilities, inlet control actuators, and geometric variations. However, it does not 
readily lend itself to test closed loop control algorithms. Requirements for inlet integration with other models 
and control algorithm verification prompted development of a nonlinear inlet model using Simulink software 
to focus on control algorithm verification and integration with other models. A more detailed description 
and validation of the Simulink inlet model is forthcoming in Kopasakis. 20 

The inlet Simulink simulation was first setup to model the same 134 grid points used in LAPIN. However, 
long simulation times motivated the need for a reduction in the grid points to reduce simulation time, while 
maintaining a normal shock position to agree with LAPIN calculations within 1.5%. This resulted in a 
model with 87 grid points, the majority of which are in the supersonic portion of the inlet. The upstream 
flow properties are static pressure, static temperature, and Mach number. An exit Mach number is used for 
the downstream boundary condition. To demonstrate the Simulink model’s ability to match LAPIN results, 
the two models were simulated at a cruise operating point. A free stream operating condition was chosen 
representing 18,288 m (60,000 ft) cruise altitude at a flight speed of Mach 2.35 and an exit Mach of 0.45 . 

The average percent dif- 
ference between the steady 
state LAPIN results and the 
Simulink simulation across all 
grid points are provided in Ta- 
ble 1 for both the 134 and 
87 grid point models. It 
can be seen that the Simulink 
simulation is comparable to 
LAPIN when using the same 
grid points and only intro- 
duces a worst- case average for 
the static pressure calculation 
of a 3% difference with re- 
duced grid points. Steady 
state plots of the static pres- 
sure, static temperature, den- 
sity, and Mach number can 
be seen in Fig. 5. The ax- 
ial location of the inlet starts 
at the “0 m” location repre- 
senting the start of the exter- 
nal compression. The external 
compression is assumed to oc- 
cur much faster than actuator 

movement, and thus will not be considered to impact the simulation dynamics. This allows only oblique 
shock relations to be used in the region upstream of the cowl lip. The CFD is used for the internal portion 
of the inlet starting downstream of the cowl lip located at “0.25 m”. The AIP is located at “1.45 m” where 





Figure 5. Inlet steady state plots of static pressure, static temperature, 
static density, and mach number for the Simulink software model compared 
to LAPIN simulation results across inlet axial stations. 
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the engine face would be located. The plots in Fig. 5 show a high level of agreement between the Simulink 
model and LAPIN across all axial locations of the inlet, with the data traces lying nearly on top of each 
other. 

Table 1. Average percent difference between LAPIN simulation results and Simulink inlet simulation results 
using 87 and 134 inlet grid points. 


Simulink 
Inlet Grid 

Static 

Pressure 

Static 

Temperature 

Static 

Density 

Shock 

Position 

87 points 

3.07 % 

0.94 % 

2.15% 

1.38% 

134 points 

0.80 % 

0.51 % 

0.58 % 

0.63 % 


A Bode plot is used to il- 
lustrate the dynamic simula- 
tion results of the LAPIN and 
Simulink simulations in Fig. 6. 

These plots show a high level 
of agreement between the two 
models. The Bode plot is 
for an exit static pressure re- 
sponse resulting from a 3% 
amplitude sinusoidal logarith- 
mic free stream pressure dis- 
turbance, chosen for illustra- 
tive purposes. The sinusoid 
starts at 20 Hz and ends at 
200Hz after ten seconds of 
simulation time. The phase 
and coherence of the two mod- 
els are almost identical, while 
the magnitudes of the Bode 
plots starts to show a small de- 
viation at frequencies greater 
than 100 Hz. 

B. Engine Model 

The engine comprises the other major component of the overall propulsion system, and is connected 
just downstream from the inlet. Initially, the modeling effort in this paper will use a turbojet engine similar 
to the General Electric J-85 for verification. More relevant supersonic engine concepts will be considered in 
the following section. The J-85 lends itself to a straightforward model development and verification because 
it has not only undergone extensive testing at NASA Glenn Research Center 24 , but many of the required 
performance parameters are accessible as well. 

The turbojet model uses a single lumped volume for each of the major components such as compressor, 
combustor, turbine, afterburner, and nozzle. The modeling approach is outlined in Seldner 25 and refined 
in Kopasakis. 3 The high level equations are defined here for completeness and understanding the level 
of fidelity of the nonlinear simulation, which is integrated into the overall propulsion model. Most of the 
geometric information for the engine is obtained from previous work done by Tesch. 26 Each of the fluid 
flow components is modeled using a set of derived conservation equations modified from the standard Euler 
form and written as continuity, momentum, and energy, Eqs. (15) to (17), respectively. These equations are 
integrated numerically using a time marching scheme and the Seldner differencing technique of the spacial 
terms. 3 


Simulink Model - - - LAPIN 

2r 

00 

■O 

.c 0— 


20 40 60 80 100 120 140 160 180 200 

Frequency - Hz 



Frequency, Hz 


Figure 6. Inlet model using Simulink software exit pressure bode plot for 
a logarithmic free stream pressure disturbance of 20 Hz to 200 Hz and the 
resulting coherence versus a LAPIN simulation bode plot. 
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A compressive component modeling 
approach is discussed in more detail in 
this section. This is done to provide con- 
text for the modifications of the bound- 
ary conditions required for inlet-engine 
model integration. Integration will be 
discussed subsequently. The three state 
variables chosen for the model of each 
component are the static density, static 
density times the total temperature, and 
mass flow rate. Since this is a sub- 
sonic system, information from two of the 
states must travel downstream, while the 
other travels upstream. Here the mass 
flow rate information travels upstream 
and the temperature and pressure travel 
downstream. A general schematic of the 
overall modeling approach for a compres- 
sive component is shown in Fig. 7. This schematic is analogous to the other major components, with the 
exception of various performance maps, which may be required to meet characteristic flow parameters. The 
variables are left generic and the subscript “n” is associated with the current volume component. The con- 
ditions of the engine component upstream are used with performance maps to generate characteristic values 
for this particular engine design. The performance map is actually two table lookups in the model to obtain 
a given pressure ratio and efficiency, which are then used to get the characteristic pressure and temperature. 
The characteristic values are then used with the associated variation of the engine governing Eqs. (15) to 
(17) to obtain the state values that are passed to the subsequent engine component. 




Figure 7. 
schematic. 


Turbine engine compressor component modeling 


Next Three Generations (N+3) Engine Design Concept 


The current engine design concept of the 
supersonic propulsion system is a Vari- 
able Cycle Turbofan Engine (VCE). The 
overall modeling approach is the exact 
same as the turbojet described above, 
with the exception that there is a more 
complicated flow path. The flow coming 
from the inlet goes through the fan com- 
ponent and is split into three gas paths 
as illustrated in Fig. 8. The primary gas 
path is nearly identical to the turbojet 
and goes through the core of the engine. 
The secondary gas path is similar to a 
typical turbofan engine in that a large 
amount of the flow is bypassed around 
the core of the engine. This secondary 
flow is mixed downstream of the turbine 
and exhausted out a common nozzle with 
the core flow. The more advance concept 


Bypass VCE 


VCE Nozzle 


Bypass - 


Customer Bleed 


Cooling Air 


Inlet 


Fan 


i r " 


i i 


► HPC — ► Burner -► HPT 

— ► LPT 

1 1 

■ HP Shaft — - 1 
1 

1 

i 

i 

i 

Parasitic Load - 1 

i 

i 

I_p Shaft 

i 


Mixer — ► Nozzle 


Figure 8. N+3 engine design concept for a variable cycle engine. 

part of this engine is a third gas path that exhausts through a 
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different nozzle. This path provides a lower exit velocity that could be used as a noise shield for the flow 
exiting the main nozzle. While Fig. 8 shows all of the major engine components that have their own conser- 
vation equations, in between each component is also an associated connecting duct that is not shown in the 
figure for simplicity. Each of these ducts are also modeled using the governing conservation equations. 

The initial conditions and geometry information to initiate the N+3 VCE comes from a design concept 
implemented in the Numerical Propulsion System Simulation (NPSS). A cruise operating point was chosen 
at 15,240 m (50,000 ft) and a Mach number of 1.6. A choked nozzle flow equation is used for the exit 
boundary condition. The average steady state percent difference between the NPSS model and the Simulink 
model is illustrated in Table 2 for the state parameters across all of the engine components. The main 
parameter of interest for this study is the thrust, calculated here as the gross thrust, which deviates from 
NPSS calculations by only 1.04%. Thrust is important because it will be the feedback parameter when the 
propulsion system is integrated with the vehicle model. 

Table 2. Average percent difference between variable cycle turbofan engine Simulink simulaton results and 
NPSS model results across all components and ducts. 


Total 

Pressure 

Total 

Temperature 

Mass 

Flow 

Gross 

Thrust 

0.97 % 

1.05% 

1 .67% 

1.04% 


C. Inlet and Engine Integration 

A complete propulsion system simulation has the nonlinear inlet model coupled with the nonlinear dynamic 
engine model. The integrated nonlinear propulsion system used in this study includes the simulated 2DB inlet 
and a representative J-85 engine models. This allows for the development of the modeling approach, control 
concepts, and interfaces for the overall propulsion system integration. Ultimately an external compression 
inlet needs to be developed and integrated with the N+3 VCE to have a more accurate picture of how 
atmospheric disturbances will impact thrust for the current supersonics design concepts. 

The integration approach is to keep the engine model input boundary condition the same, thus expecting 
a total pressure and total temperature. The AIP boundary condition of the inlet is removed and replaced 
with a very small duct volume that can be integrated into the engine simulation. The duct volume is modeled 
similarly to the ducts in the N+3 VCE engine, where the conservation equations become Eqs. (18) to (20). 


dt 


(psT t ) 


1 _ 

V 


j t {p s ) = - m en9ine ) 

(18) 

d A ■ 

(777,) — p'i r nlet ex it ^ 

(19) 

(rpinlet exit- 1 rfrfnlet exit _ pinlet exit ^engine^ 

(20) 


The volume that includes the AIP uses its calculated states to output the total pressure and total temperature 
required for the engine. In addition, the removal of the Mach number AIP boundary of the inlet requires 
information of the engine mass flow rate to be accepted as the engine feedback information using standard 
gas property relations. 

The superscripts used in the above equations indicate the location of the variable relative to the AIP and 
engine face. These equations are then tied together using the state equation, Eq. (12), to obtain the inlet 
exit state variables for the downstream engine face and the upstream inlet grid points. 

D. Disturbance Models 

The primary source of disturbances investigated in this study are upstream flow perturbations due to 
atmospheric wind gusts, using the atmospheric model outlined in Kopasakis. 11, 12 An advantage of this 
atmospheric model is that it allows for time domain simulation of disturbances that are fractional order. 
These disturbances are used in determining the effectiveness of the control to limit thrust oscillations and, 
in the case of the mixed compression inlet, avoid an unstart condition. For this control design, an important 
goal is to assess the worst possible disturbance that the propulsion system is expected to encounter during 
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operation. Due to the difficult nature of defining a worst-case for atmospheric disturbances, a representative 
severe-case scenario is obtained from an approach where several sinusoidal wind gust waves are summed to 
obtain a repeating pattern that can be used to test the control algorithm. The highest frequency wind gusts 
based on the vehicle speed and altitude are on the order of tens of hertz. 

E. Controller 

The engine control design used in this study has aspects of classical loop shaping 13 design and Quantita- 
tive Feedback Theory. 27 The engine controller for the shaft speed was first designed on linear models, 15 and 
later implemented in the nonlinear model turbojet and the N+3 VCE model with a slight reduction in the 
gain to improve stability. The shaft speed is controlled using the engine fuel injector that has a bandwidth 
of 6 Hz. In addition, the turbojet and VCE engine design both have variable nozzle and compressor guide 
vanes that are set on schedules and not actively controlled. The bandwidth of approximately 6 Hz for a fuel 
injector will not be able to attenuate disturbances in the tens of hertz. This results in the inlet bypass door 
being the primary means to attenuate the higher frequency disturbances and reduce thrust oscillations 28 ’. 29 

A controller for the 2DB inlet normal shock position was developed using the same approach based on 
linear models obtained from LAPIN. 28 The control design uses a bypass door actuator with a bandwidth 
of approximately 175 Hz, which is modeled using a simple transfer function. 15 However, the control design 
from the previous work encountered stability issues when implemented in the Simulink model. This could 
be a result from the change in the model’s operating point combined with the change to control the normal 
shock closer to the throat and thus altered the linear operating point from which the controller was designed. 
Therefore an alternate PID controller with a low pass filter was designed, but with a lower bandwidth. The 
transfer function of the PID controller is listed in Eq. (21). 

1 165 5 

C in i e t = 0.0148 + 6.0482- - (3.7 x 1CT 5 )- — — r (21) 

s 1 + 165.5^ 

The controller’s lower bandwidth results in an inability to attenuate high frequency disturbances (on the 
order of tens of hertz) that could be encountered in either the atmospheric or ASE disturbances. While this 
is a concern, the emphasis here is propulsion system integration to provide a platform for control design, 
rather than a detailed control development. Also, current concepts for use of an external compression inlet 
will not have a bypass door due to the lack of a need to control the normal shock and thus will not be able 
to actively attenuate the high frequency disturbances as well. 

IV. Results 

The purpose of these simulation results is to demonstrate the integrated propulsion system model with 
an included controller. The example simulation results investigate the impact of atmospheric disturbances 
on the normal shock position and thrust oscillations from the propulsion models about a cruise operating 
point. Presentation of the results is divided into three sections that represent the main propulsion system 
platforms. The first section will investigate the 2DB inlet model alone for its ability to use the bypass door 
to control the normal shock position in the inlet, second the integrated 2DB inlet and representative J-85 
engine model will be used to investigate a fully integrated supersonic propulsion system’s thrust oscillations, 
and finally a preliminary look at a current concept VCE engine model will be used to investigate thrust 
oscillations. 

A. Inlet Results 

The Simulink based 2DB inlet model is simulated at a cruise condition at 18,288 m (60,000 ft) and a flight 
Mach number of 2.35. The AIP boundary condition is set to a Mach number of 0.45. The normal shock is 
commanded at the 1 m location of the inlet for illustrative purposes. To provide a visual reference Fig. 2 
can be used where 0 m would refer to the start of the inlet, 0.25 m would refer to the start of the internal 
section of the inlet, 0.66 m refers to the minimum area or throat of the inlet, and the AIP would be at 1.45 
m. The closed loop controller moved the normal shock position upstream and downstream from its nominal 
1 m location due to a commanded position change of db 3% from nominal, illustrated in Fig. 9. The normal 
shock position is shown to have a very good response with no overshoot. While the ability to command 
the shock to the desired location provides some confidence in the PID controller, the goal is to attenuate 
atmospheric disturbances. 
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The atmospheric turbulence model is used at the 
cruise condition to provide disturbances to the input 
boundary conditions of temperature, Mach number, 
and pressure. The boundary condition disturbances 
are illustrated in Fig. 10. A sum of sinusoids cre- 
ates repeating disturbances that are representative 
of the one-second time slice shown in Fig. 10. The 
main thing to take away from this plot is the large 
drop in boundary condition values at the six second 
mark. This atmospheric disturbance was used for 
both the inlet simulation results and the integrated 
inlet-engine simulation results. 

The atmospheric disturbance results in a large 
normal shock movement for the open loop inlet as 
seen in Fig. 11. The open loop normal shock re- 



Time, s 

Figure 9. Inlet normal shock position command step 
changes controlled by inlet bypass door actuation. 


sponse oscillates about the 1 m location; however, the large drop in free stream conditions at the six second 
mark causes the normal shock to move to the throat location of 0.66 m. This is a marginally stable location 
for the normal shock, meaning that even slightly larger perturbations could cause the normal shock to be 
ejected from the inlet, causing an unstart condition. This is a dramatic event that results in sudden loss of 
thrust, and could create large ASE modes. 



Time, s 


Figure 10. Simulated free stream perturbations based on atmospheric disturbance model wind gusts 

The introduction of the closed loop normal shock 


control greatly attenuates the affects of the atmo- 
spheric disturbance as shown in Fig. 11. The con- 
trolled bypass door attenuates the normal shock 
movement throughout the one-second repeating pe- 
riod of the atmospheric disturbances. The large 
transient that nearly caused the open loop response 
to unstart only causes the controlled normal shock 
response to move 0.08 m upstream from the com- 
manded position and 0.1 m in the relatively benign 
downstream direction. This control design for the 
bypass door, provided a very desirable response to 
the disturbance. Further development of the con- 
troller could allow for the normal shock to be com- 
manded closer to the throat instead of the chosen 
1 m illustrative location. Also, given that these are 
atmospheric disturbances, it is difficult to determine 
exactly the worst-case disturbance, thus the disturbances here are best classified as severe. 



Figure 11. Simulated inlet open-loop and closed-loop 
normal shock position response to an atmospheric dis- 
turbance. 
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B. Integrated Inlet and Turbojet Engine Results 


The integrated model investigates the closed loop response of the 2DB inlet simulation and representative 
J-85 engine model to the atmospheric disturbance described in the inlet section. The operating condition is 
representative of a cruise condition at 19,812 m (65,000 ft). This condition was chosen for the operating point 
as it represents a wind tunnel test condition in which the specified components were tested in the NASA 
Glenn Research Center 10x10 supersonic wind tunnel. The integrated model only focuses on the closed loop 
response, because the open loop response that caused the inlet to nearly unstart introduced oscillations that 
can cause the integrated model to have instabilities. The inlet normal shock position is still commanded 
to the lm location as can be seen in Fig. 12. The closed loop response of the integrated model has an 
advantage of having both a bypass door and fuel injector actuation for control purposes. The advantage can 
be seen as the normal shock position movement is reduced in the integrated model over that of just the inlet 
model alone. The integrated closed loop response has a maximum deviation of 0.01 m and the inlet model 
had a maximum deviation of 0.1 m, approximately. This order of magnitude better attenuation illustrates 
a more rigorous control design coupling the inlet and engine control algorithms could provide a substantial 
improvement. 




Figure 12. Integrated inlet and engine simulation, Fi S ure 13 ‘ Simulated propulsion system net thrust os- 
with Closed-loop normal shock position control re- cillations due to atmospheric disturbances with closed- 
sponse to atmospheric disturbance. loo P normal shock P osition and rotor s P eed control. 


Finally, the integrated model is used to show the closed loop net thrust oscillations due to an atmospheric 
disturbance, illustrated in Fig. 13. The nominal resulting net thrust from the integrated model is 6,900 N 
(1551 lbf). The larger sharp decrease in the free stream conditions that occurs at six seconds causes a large 
thrust oscillation. The mach number disturbance is plotted on the secondary y-axis to provide a reference 
of the time delay from the input disturbance to the resulting thrust oscillation. The larger transients cause 
an oscillation of approximately 1,100 N ( 247 lbf) or a 16% sudden change in the thrust of the propulsion 
system. This is a significant perturbation that warrants further investigation into the propulsion system 
effects on the vehicle. 

C. N+3 Variable Cycle Turbofan Engine Results 

The integrated model described in the previous section provides valuable modeling methodology verification 
and a preliminary test bed; however, the goal is to have an integrated model for the current design concepts 
of the supersonics project. In this section the N+3 VCE engine results will be presented to test out the 
controller and examine a severe-case scenario for atmospheric disturbances with a potential to cause thrust 
oscillations at a cruise operating point. The engine controller component provides a reasonable initial look at 
the thrust oscillation suppression as opposed to including an inlet control, because the inlet being suggested 
for the supersonics project is an external compression inlet and will not have any of the inlet control actuators 
previously described. The external compression inlet will provide a large damping volume for the disturbances 
to attenuate before reaching the engine face, but no active control mechanisms are currently planned to help 
attenuate the disturbances. Thus, in this study the use of the engine model alone will provide a worst-case 
configuration, in the sense that without the inlet no damping volume is used to attenuate the atmospheric 
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disturbances prior to reaching the engine face. Only the engine fuel injector that controls the fan speed will 
be used to dampen out any oscillations. 

Results from the N+3 VCE are all ob- 
tained at the operating point of 15,240 m 
(50,000 ft) and Mach 1.6. The control 
performance was first examined by com- 
manding the fan speed through a step up 
and step down transient with a 3% am- 
plitude change. The results of this first 
test is shown in Fig. 14. The closed loop 
fan speed response illustrated in the up- 
per plot of Fig. 14 shows an ability to 
track the commanded fan speed changes 
well with only small underdamped oscil- 
lations before settling on the steady state 
value The lower illustration in Fig 14 Figure 14. VCE fan speed command step changes controlled by 
i i -a . . r i fuel injector actuation and the resulting gross thrust change. 

shows that a 3% change m the fan speed 

results in a 3.2% change in the gross thrust produced by the engine from a steady state value of 44,100 N 
(9,914 lbf). 

The atmospheric turbulence model is applied to provide disturbances to the input boundary conditions of 
total pressure and temperature at the engine face. This is very similar to the boundary condition disturbances 
that were applied to the inlet, which are illustrated in Fig. 10. The distinction is that the velocity component 
of the disturbance that was used to excite a Mach number disturbance is not used for the input disturbance 
to the engine. The engine model only has the two upstream boundary conditions whereas the inlet has three 
upstream boundary conditions pertaining to the disturbances in Fig. 10. The Mach number disturbance 
could have been converted into an additional temperature disturbance, but was not considered for this 
study. Similar to the inlet atmospheric disturbance, a sum of sinusoids creates a repeating pattern that is 
representative of the one-second time slice. Again, the disturbance magnitude at the 6.0 second mark is a 
significant disturbance to the propulsion system. 
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Figure 15. N+3 closed-loop gross thrust response to atmospheric 
disturbances at engine face, with a fan speed controller. 


The engine face disturbance results 
in a thrust oscillation about the nominal 
44,100 N (9,914 lbf) value of only about 
2.76% as shown in Fig. 15 prior to the 
larger transient at six seconds. During 
the transient at six seconds where there 
is a larger drop in the free stream vari- 
ables, an overall gross thrust oscillation 
of 3,980 N (895 lbf) or a variation of 9% 
occurs. This is over a short period of ap- 
proximately 0.1 s. The pressure distur- 
bance is plotted on the secondary y-axis 
to provide a reference of the time delay 
from the input disturbance to the result- 
ing thrust oscillation. Thrust variations 
experienced in this study are expected to 
be reduced once the external compression 


inlet is integrated with the engine; however, this preliminary investigation does suggest that the thrust os- 
cillation is significant enough to merit further study. 


V. Conclusion 

An integrated inlet and engine propulsion system component model has been developed and integrated 
with a control algorithm, where the inlet and engine components have previously been developed. The 
integrated propulsion system model presented in this work is suitable for incorporation into an overall 
supersonic vehicle aero-propulso-servo-elastic (APSE) model. The simulation includes an inlet normal shock 
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position control and an engine rotor speed control for the propulsion system of a supersonic commercial 
transport. An atmospheric turbulence model is used to perturb the integrated propulsion model to illustrate 
system stability and expected thrust oscillations. The integrated propulsion system results illustrate that 
the inlet and engine controller allow for stable operation without a threat of an unstart condition for the 
single operating point investigated in this study. However, the thrust variation experienced has the potential 
to excite vehicle aero-elastic modes. The thrust oscillation obtained from the atmospheric disturbances was 
approximately a 15% variation from nominal. Such an oscillation will cause the aero-elastic modes of the 
vehicle to be impacted and further study of the coupling between these dynamics is required. A preliminary 
investigation of thrust oscillations induced by a severe atmospheric disturbance using a current supersonics 
project propulsion system concept, the next three generations (N+3) variable cycle engine (VCE), without 
an inlet model was conducted. The thrust oscillations varied by approximately 10% about the nominal value, 
due to the disturbance. These results are expected to be a worst-case configuration, because the disturbance 
is placed directly at the engine face; however, there is potential for more severe atmospheric disturbances 
to be applied. The 10% magnitude thrust oscillation from the N+3 VCE design concept could excite the 
aero-elastic modes of the vehicle and further study is required. 
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